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Abstract 

In this paper, the branching ratios and the direct CP-violating asymmetries for decays — t- 
/o(980)^**^ and B~ — /o(980)A'*^ by employing the perturbative QCD factorization approach 
are studied. In the two-quark model supposition, /o(980) is commonly viewed as a mixture of ss 
and nn{= {uu + dd) / , that is |/o(980)) = \ss) cos 9 +\nn) sin 6, where 6 is the fo — cr mixing 
angle. We find that the non-factorizable /o emission type diagrams can give large contributions 
to the final results, which are consistent with the present experimental data and the upper limit 
in the allowed mixing angle ranges. We predict that the direct CP asymmetry A'(^p{fo{980)K*^) 
is small, only a few percent, which can be tested by future B factory experiments. 
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I. INTRODUCTION 



In order to uncover the mysterious structure of the scalar meson /o(980), intensive 
studies have been done since it was firstly observed in vrvr scattering experiments [H- 
There is still no consensus on the essential inner structure of /o(980). Some people 
consider it as qq state [2I], or qqqq four-quark state jl], other people think that it is not 
made of one simple component but might have a more complex nature such as having a 
KK component 0, 0], or mixing with glueball [6l-[8|, or even superpositions of the two- 
and four-quark states. 

The B decays involved in the /o(980) in the final states are studied by employing various 
factorization approaches, such as the generalization approach [9], the QCD factorization 
(QCDF) approach [lol- the perturbative QCD (PQCD) approach |13l- 16|. In these 



calculation, the scalar meson is usually viewed as a mixture of ss and nn{= (uu+dd) / a/2), 
that is 

|/o(980)) = |ss)cos^+ |nn)sin^, (1) 
where 6 is the fo — (t mixing angle. About the value of 6, there are many discussions in 



the phenomenal and experimental analyses [17|, |18|. But unfortunately it is difficult to 
find a unique mixing angle to describe /o — o" mixing. 

On the experimental side, for /o(980) emerging as a pole of the amplitude in the S 
wave [l9[, many channels such as B ^ /o(980)-ft' can be obtained by fitting of Dalitz 
plots of the decays B — t- -n-^n^K and B — )■ KKK and so on |20-24|. Although many 



such decay channels that involved /o(980) in the final states have been measured over the 
years, it has yet not been possible to account for the this scalar meson inner structure , 
i.e. whether one deals with a two- or rather a four-quark composite, because there still 



lack precise and enough data. For our considered decays, the measured values are 25| : 

Br{B- fo{980)K*-) = (10.4 ± 2.6) x 10-^ (2) 
Br{B^ /o(980)i?*°) < 8.6 x 10"^ (3) 

It is noticed that we have assumed i?r(/o(980) — ?■ tt+tt^) = 0.50 to obtain the upper 
experimental branching ratios. 

In this paper, we will study the branching ratios and the direct CP asymmetries of 
B° /o(980)^*° and B~ /o(980)i^*" within perturbative QCD approach based on kr 
factorization. In the following, /o(980) is denoted as /o in some places for convenience. It 
is organized as follows. In SectHIl the relevant decay constants and light-cone distribution 
amplitudes of 5, /o and i^'*are discussed. In Sec JIIIl we then analysis these decay channels 
using the pQCD approach. The numerical results and the discussions are given in section 
IIV[ The conclusions are presented in the final part. 

II. DECAY CONSTANTS AND DISTRIBUTION AMPLITUDES 

Now we present the wave functions to be used in the integration. For the wave function 
of the heavy B meson, we take: 

b) = + mBh5(pB{x, h). (4) 



2 



Here only the contribution of the Lorentz structure (j)B{x,b) is taken into account, since 
the contribution of the second Lorentz structure (pB is numerically small [26| and has been 
neglected. For the distribution amplitude (pB^x^h) in Eq.(jl]), we adopt the model 

M^M = Nbx\1 - exp[-^ - (5) 

where Ub is a free parameter, and the value of the normalization factor is taken as A^^ = 
91.745 for Uf, = 0.4 in numerical calculations. 

In the two-quark model, the vector decay constant ff^ and the scalar decay constant 
//q for scalar meson /o can be defined as: 

(/o(p) 1927^91 1 0) = ffoPt,, (6) 



{fo{p)\q2qi\0) = nifjf^. (7) 

Owing to charge conjugation invariance or the G parity conservation, the neutral scalar 
meson /q cannot be produced via the vector current, so ff^ = 0. Taking the mixing into 
account, Eq.([7]) is changed to 

{fS\dd\0) = {fS\uu\0) = -^rnjJl, {f^\-ss\0) = mj^. (8) 

Using the QCD sum-rule method, one can find that the scale-dependent scalar decay 
constants f^^ and /^^ are very close |ll|. So f^^ = fj^ is assumed and we denote them as 
ffQ in the following. 

The light-cone distribution amplitudes (LCDAs) for the scalar meson /o can be written 

as: 

{Mp)\qi{zU2mO) = -^J\xe-^- 

x{^^f,{x) + mf^^%{x) + mfy{yi+4^ - l)<l>J^(x)},,, (9) 

here n+ and n_ are light-like vectors: = (1, 0, Ot), = (0,1, Ot), and n+ is parallel 
with the moving direction of the scalar meson /q. The normalization can be related to 
the decay constants: 



dx^f,{x)= I rfx$J„(x) = 0, / dx^%{x) = -i^. (10) 
The twist-2 LCDA can be expanded in the Gegenbauer polynomials: 

^ oo 

$^„(a;,/x) = ^-^4(/x)6x(l-a;)5^i?„(/i)C^/2(2a;-l), (11) 



m=l 



the values for Gegenbauer moments -B3 have been calculated in [11 



as: 



5i = -0.78 ±0.08, 53 = 0.02 ±0.07. (12) 
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These values are taken at /i = 1 GeV and the even Gegenbauer moments vanish. 

As for the twist-3 distribution amphtudes 
the hterature, so we adopt the asymptotic form : 



As for the twist-3 distribution amphtudes and they have not been studied in 



1 



1 



For our considered decays, the vector meson K* is longitudinally polarized. The lon- 
gitudinal polarized component of the wave function is given as: 



1 



(14) 



where the first term is the leading twist wave function (twist-2), while the second and 
third term are sub-leading twist (twist-3) wave functions. They can be parameterized as: 



K* 



6x(l - x) 1 + ai/^.Cf/^(2x - 1) + a2K*C^r^\^x - 1 



3/2/ 



(15) 



$^.(x) 



:i-2a;), $^.(x) 



— ^^f2a;- 
2v^ ' 



(16) 



where the Gegenbauer moments a\K* = 0.03,a2_ft:* = 0.11 j27| and the Gegenbauer poly- 
nomials C^{t) are given as: 



(t) = 3t, cY\t) = ^{5t'-l), 
Cl'\t) = lt{7t'-3). 



(17) 
(18) 



III. THE PERTURBATIVE QCD CALCULATION 

Under the two-quark model for the scalar meson /o supposition, we would like to 
use pQCD approach to study B decays into /o and K*. The decay amplitude can be 
conceptually written as the convolution, 

A{B^foK*) ~ Jd%d%d% TT[CmB{ki)^f„{k2)^K<h)H{h,k2,h,t)], (19) 

where /cj's are momenta of anti-quarks included in each mesons, and Tr denotes the trace 
over Dirac and color indices. C{t) is the Wilson coefficient, which results from the radiative 
corrections at short distance. In the above convolution, C{t) includes the harder dynamics 
at larger scale than Mb scale and describes the evolution of local four-Fermi operators 

from mw (the W boson mass) down to t ~ 0{^/AMb) scale, where A = Mb — rrih- 
The function H{ki, k2, k3,t) describes the four-quark operator and the spectator quark 

connected by a hard gluon whose is in the order of AM^, and includes the 0{^/AMb) 
hard dynamics. Therefore, this hard part H can be perturbatively calculated. 
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Since the b quark is rather heavy we consider the B meson at rest for simphcity. It is 
convenient to use hght-cone coordinates {p'^,p~ , p^^) to describe the meson's momenta by 



= -^{P^ ±P^), and pt = 



(20) 



Using these coordinates the B meson and the two final state meson momenta can be 
written as 



Pi 



rriB , 



B 



1,1, Ot), P2 



rriB 



' ' ^ V2 

respectively. Here we have the mass ratios 



(l-r^,r2 Ot), B 



TUB, 2 
1 r 

V2 



rfo = mfJuiB. 



(21) 



(22) 



Putting the anti-quark momenta in B, /g, K* mesons as ki, k2, and k^, respectively, we 
can choose 



ki = (siP+, 0, kiT), k2 = ix2P^, 0, ksr), ks = (0, 0:3^3", ksT 



(23) 



For these considered decay channels, the integration over , 5 ciiici in eq.( fT9l) will 
lead to 

A{B foK*) ~ Jdxidx2dxsbidbib2db2b3db3 

■Tr [C(t)$B(xi, 6i)$/o(x2, &2)$x*(x3, bs)H{x„ h, t)St{x^) e-^(*)] , (24) 

where bi is the conjugate space coordinate of kix, and t is the largest energy scale in 
function H{xi,bi,t). In order to smear the end-point singularity on Xi, the jet function 
St{x) [2^, which comes from the resummation of the double logarithms In^Xj, is used. 
The last term e~'^^*^ in Eq. (l2^ is the Sudakov form factor, which suppresses the soft 
dynamics effectively (29| . 

For the considered decays, the related weak effective Hamiltonian H^fj can be written 
as [3C 



n 



eff 



71 



q=u,c 
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1=3 



(25) 



with the Fermi constant Gp^ = 1.16639 x 10 ^GeV ^, and the CKM matrix elements V. 
We specify below the operators in T-Leff ^01 b ^ s transition: 



(26) 



where a and /3 are the S'[/(3) color indices; L and R are the left- and right-handed 
projection operators with L = (1 — 75), R = (1 -I-75). The sum over q' runs over the quark 
fields that are active at the scale /i = 0{mb), i.e., {q'e{u, d, s, c, b}). 



0^ 




Up-f^Lba , 


0^ 


= SaYLUa 


■ up'^^Lbp , 

Eg' Q.'plM'a , 
Eg' Q'pll^RQ'a , 


03 


= Sal^Lba ■ 


Eg' Q'lS^i^Lq'f^ , 

Eg' Q'/Slt^Wp , 


04 


= SaYLb,3 


05 


= Sa'y^Lba ■ 




= SaYLbfS 


07 




■ Eg' (iq'Q'pli^RQ'p 

■ Eg' eqrq'.^M, 


Os 


= ^SaYLb^ 


■ Eg' ^q'Q'fS^t^RQ'a 

■ Eg' (iq'Q'plM'a 


09 


= ^SaYLba 


Oio 


= \sa-i''Lbf, 
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FIG. 1: Diagrams contributing to the decay — )• fo(980)K 
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In the following, we take the — )■ foK*^ decay channel as an example to expound. 
There are 8 type diagrams contributing to this decay, as illustrated in FigJH For the 
factorizable emission diagrams (a) and (b), operators Oi_4^9^io are {V —A){y —A) currents, 
and the operators Os.g have the structure of {V — A){y + A), the sum of the their 
amplitudes are written as FeK* and i^e^*, respectively. For (/o|(77/x9|0) = 0, one then 
finds that 



eK* 



0. 



(27) 



In order to get the right flavor and color structure for factorization to work, a Fierz 
transformation for the {V — A){V + A) operators may sometimes be needed and then the 
corresponding amplitude is 



77.P2 



pi roc 

Jo Jo 



xEei{t)KiXi, X3, 61, 63) - 2rK*^K*ix3)Eeiit')K{x3, Xi, 63, bi] 



(2^ 



where Cp = 4/3 is the group factor of the SU{3)c gauge group. The functions Eei{t^'^ 
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and hp including the Sudakov factor and jet function have the same definition as those in 



Ref. 15, 16 



For the non-factorizable diagrams 1(c) and 1(d), all three meson wave functions are 
involved. The integration of 63 can be performed using the S function 5(63 — 62), leaving 
only integration of 61 and 62- For the {V - A){V - A), {V -A){V + A), {S - P){S + P) 
operators, the results are 

■MeK* = 327rCirm^/ \/2Nc / dxidx2dxs / bidbib2db2^B{xi,bi)^ f^^(x2) 

Jo Jo 

[{X2 - l)$i^.(x3) - r^.a;3($^.(x3) + (x3))]E^(^)/l„(xl, 1 - X2, X3, bi, 62) 
-[(X2 +X3)$i^.(x3) +r^.X3 ($^.(X3) -$^.(X3)] E'^i{t')K{xi, X2, x^MM) [,(29) 



MeK" — ^'2'^Cpm%r f^/ \/2Nc / dxidx2dxs / bidbib2db2^B{xi,bi) 

Jo Jo 

x|£;^(t)/^,„(Xl,l -X2,X3,6i,62) [{x2 - 1)^k*{x3) {^i{^2) - *Jo(^2)) 

-rK4x2 - 1) i^U^s) + ^U^s)) {^iix2) - ^lix2)) 
+rK-^X3 {<t>Ux3) - <^UX3)) {<^1{X2) + $^0(^2))] 
+E'^i{t')hn{xi,X2,X3,bi,b2) [X2^K'{X3) (^/^ (^2) + (^2)) 
-rK'X2 ($^.(X3) + ^K*ixs)) ($f„(x2) + ^/^(xa)) 

-TK^Xs i^Uxs) - ^UX3)) {^■UX2) - $J„(X2))] ) , (30) 



/* 1 POO 

M-eK* = -^2-nCFm%/ ^2Nc / dxidx2dx3 / 61 62*^^2 $B(a;i, fei)$/o (0:2) 

Jo Jo 

x| [(-a;2 + 3:3 + 1)^k4x3) + rK*X3 {^K^ixs) - ^K'ixs))] 

XE'^i(t)hn{Xi, 1 - X2, X3, 61, 62) - E'^i(t')hn{Xi,X2, X3, 61, 62) 

X [X2<^K'{X3) + rK'X3{^K*ix3) + '^K*{X3))] | • (31) 

For the non-factorizable annihilation diagrams (e) and (f), again all three wave 
functions are involved. MaK*, ^ak* ^^"^ ^aK* describe the contributions from the 
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{V -A){V -A), {V - A){V + A) and {S - P){S + P) type operators, respectively, 

f'\ f'OC 

MaK* = 327rCj.myy2iVc / dxidx2dx^ I hidhih2dh2^B{xi,hi) {[x2^K*{x^)^f^{x2) 

Jo Jo 

+rK*r/o$J^(x2) {{X2 + X3 - 1)$^*(X3) + {X2 -X3- 1)$^*(X3)) 
-rK^rf,^%(x2) {(X2 -xs + S)^%.(xs) + (X2 + Xs - 1)$*^.(X3))] E'Jt) 

X2,Xs,bi,b2) + E'^-{t')h'^^{xi,X2,X3,bi,b2) [{xs - l)^K*{x3)^fo{x2) 
+rx*r/o$^^(x2) {{X2 -X3 + 1)^%.{X3) - {X2+X3 - 1)$^*(X3)) 

-rK*rf,^l{x2){{l-X2-Xs)^'K*ixs) + {l + X2-Xs)¥j,.{x2))]} , (32) 

Mak* = ^'27rCFm%/y/2Nc / dxidx2dx^ I bidbib2db2^B{xi,bi) 

Jo Jo 

x| [rK^il + Xs)<^>f,{x2){<^>'K^{Xs) + $^.(0:3)) - rf,{x2 - 2)<^>K*{X3) {<^%{X2) 

~*/o(^2))] E'^i(t)hna(xi,X2,X3,bi,b2) - E'^i{t')h'^^{xi, X2, x^MM) Vk* 
x(x3 - l)$;„(a;2)(<f^.(a;3) + ^'k*{x,)) - rf,X2^K*{x,){^%{x2) - ^l{x2))\ |,(33) 

dxidx2dx3 / bidbib2db2^B{xi,bi) 
Jo Jo 

x| [(x3 - 1)^/0(3:^2) ^K*ix3) + 4:rK*rf^^%{x2)^K*ixs) + rK*rfo {{x2 - X3 - 1) 

X ($^.(X3)<I'J„(X2) - <^'k,{x3)<^1{x2)) +(x2 +X3 " 1) ($^. (xa)^ J, (^2) 
-$*;,^,(,T3)$J,(X2))] ^^i(t)/l„a(2;i, Xs, Xa, 61, 62) + E'^i{t')h'^^{xi,X2, X3, 61, 62) 
X [x2^fo{x2)^K*{x3) - X2rK*rf,{^%{x2) - (^2)) (*^. (X3) + ^^^(xa)) 

-r^.r/„(l - X3)(^%(X2) + ^l(x2))(^U^s) - ^U^^))] } ■ (34) 

The factorizable annihilation diagrams (g) and (h) involve only two final state mesons' 
wave functions. There are also three kinds of decay amplitudes for these two diagrams. 
FaK* is for (y — A)(y — A) type operators, F^^» is for {V — A)(y + A) type operators, 
while is for (5* — P){S + P) type operators: 

/■I poo 

FaK* = = -87rCj.m^/s / ^^2^X3 / 62(^^2 &3C^&3 {[(2^3 - l)^if*(2^3)^/o(2^2) 

Jo Jo 

-2rK*rf,{xs - 2)$^.(x3)$J^,(x2) - 2r^.r/,a;3$*x.(a;3)$J^(x2)] 

XEai{t)ha{x2, 1 - X3, 62, ^s) + Eai(t')ha{l - X3, ^2, 63, 62) 

x[x2$i,*(a;3)$/o(a;2) - 2rK*rf,<i>U^s){{x2 + l)$?„(a;2) - {x2 - l)$jj]} ,(35) 
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/• 1 noo 

Jo Jo 



x|[ri,.(x3 - 1)<I'/0(X2)($^*(X3) -$^,(X3)) + 2ry„<l>K,(x3)$j3(x2)] 

xEai(t)ha{x2, 1 - ^3, 62, ^s) " Eai{t')ha{l - X^,X2MM) 

X [2r;,. <|.^. (X3)<f /o (X2) - r^oa^2$i.* {x,) ($ J„ (X2) + $f„ (X2))] } • (36) 

If we exchange the K* and /o in FigJH the corresponding expressions of amphtudes for 
new diagrams will be similar with those as given in Eqs.f l29tl36|) and can be obtained by 
the replacements: 

*/o ^ (x) , ^% (x) M (x) , (x) M (x) , r/„ M r;^* , (37) 

since the wave functions for the mesons /o(980) and K* have exactly the same form. 
The only difference is some normalization constants for the different twist distribution 
amplitudes. That is, the factorization formulae for (a) and (b) in the new diagrams 
amplitudes are written as: 

Fef^ = —SnCF'm'^fK* / dxidx2 / bidbib2db2^B{xi,bi) 

Jo Jo 

[(l+X2)$/o(x2) -r/o(l -2X2) ($/o(x2) -$]J^(X2))] Eei{t)he{Xi,X2,bi,b2) 

-2rf,<l>%{x2)E,,{t')K{x2,x,,b2,bi)\ , (38) 



X 



= 0. (39) 

/o 



(S T) 

Since we have chosen the momentum fraction at the anti-quark, we should use $^ ' ^ 



x) and $^*'*(1— x) for the mesons /o(980) and K* in the calculation. But for simplicity, we 
use $j^'"^'*(x) and $^'*^(x) to denote $^^'^■'(1 — x) and $^f''(l — x) in the upper formulae. 

Combining the contributions from different diagrams, the total decay amplitudes for 
the decays 5° — foK*'^ and — )■ foK*~ can be written as: 

MifoK') = iuM,K'C2F,{e) - 6 {Fe^o(a4 - ^)Fi(^^) + [F^F^ + Ff^lF2(^)](a6 - |) 

+ [Me/o(C3 - ^) + MeK*(2C4 + + M^K* (C3 ^C,-^-^) 



xF2{e) + [MgF,{e) + MS.F2(^)](C5 - ^) + m5,2. 



(2C6 + y)Fi(^) 



+ [M^s^MO) + M,K*F2{e)] [C, - ^) + [M^lF,{e) 



+ iCe~Y^F2{9) 

+M^K^F2{e)] {C, - ^) + [F^f,F,{0) + FaK*F2{e)]{a, - ^) 

+ [Fa7oM0) + FriF2{e)]{a, - |)} , (40) 



TABLE L Liput parameters used in the numerical calculation 
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25| 



Masses 


nif^^ = 0.980 GeV, 


rriK* 


= 0.892 GeV, 




Mb = 5.28 GeV, 






Decay constants 


/b = 0.19 GeV, 


ffo 


= 0.37 GeV, 




fK* = 0.217 GeV, 


Ik* 


= 0.185 GeV, 


Lifetimes 


tb± = 1.638 X 10-12 


tro = 


1.530 X 10-12 


CKM 


Vtb = 1.0, 


Vts 


= -0.0387, 




Vus = 0.2255, 


Vub = 


0.00393e-*60° . 



xCi + FaK-a^) F^m - it [FefM + aio)i^i(^) + FfjlFi{9){a, + as) 



+M,MC, + C4 - Y - ^)F2{e) + MgF,{e){C, + C,) + Mfi. 



+ [M,f,F,{e) + M,K*F,{e)]{C, + C9) + [Mf^;Fi(^) + M^^,F,{e)] 
x(C5 + C7) + [FahF^iO) + F,^.F2(^)](a4 + aio) 

+ [^^5^1 W + i^ri.i^2W]K + as)} , (41) 

where = KifeKTs^^t = ^fe^s -Fi(6') = sin 6'/-\/2, ^2(6') = cos 6'. The combinations of 
the Wilson coefficients are defined as usual 13 ll: 



aM = CM + aM = CM + 

ai(/i) = Ci(/i) + ^^^^^1^, i = 3,5,7, 9, 

a,(/i) = CM) + ^ = 4, 6, 8, 10. (42) 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

In the numerical calculation, we will use the input parameters as listed in Table [B 
From Eq. fl38|) . we can find the numerical values of the corresponding form factor 



pB^^^HM)^^ maximal recoiling: 



which is smaller than _F'^'^-^/o(980)(drf) _ g^^-j? jjj^ using different values for the threshold 
parameters c in the jet function. 
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In the B-rest frame, the decay rate of 5 — )■ /o(980)-ft'* can be written as: 



where Tf^^rx* have been defined in Eg. (122!) and Jvl is the total decay amphtude of i? — )■ 
/o(980)if*, which has been given in section UTTl 

Using the wave functions as specified in the previous section and the input parameters 
hsted in Table [U it is straightforward to calculate the CP-averaged branching ratios for 
the considered decays. 

If /o(980) is purely composed of ss, the branching ratios of i? — t- /o(980)i^'* are: 

5r(5° ^ /o(980)i?*°) = X 10-^ (45) 

Br{B- ^ Um)K*-) = (IbAtlitlttti) X 10-^ (46) 

where the uncertainties are from the decay constant of /o(980), the Gegenbauer moments 
Bi and B^. If /o(980) is purely composed of nn, the branching ratios for B /o(980)ii'* 
are: 

BriB' ^ /o(980)ir*0) ^ (2.2t0:310jtJ:2) x 10-^ (47) 
BriB- ^ fo{980)K*-) = (4.0+°^;°:^+};?) x 10-^ (48) 

where the uncertainties are from the same quantities as above. 

The branching ratio for decay B^ — t- /o(980)-ft'*~ in the upper extreme case is consistent 
with QCDF results |l2|: 

I 6.9 X 10 °,for/o(980) = nn. 

The Branching ratio of i? — ?■ /o(980)i^* depends on the mixing angle 6 of strange and 
nonstrange components of the /o(980). In Fig 13], we plot the branching ratios as functions 
of the mixing angle 6. Using the above mentioned range of the mixing angle, we obtain: 

BriB- ^ fo{980)K*-) = (11.7 ~ 14.6) x 10~^ (50) 
5r(5° ^ /o(980)i?*°) = (11.2 ~ 13.7) x 10-^ (51) 

for 25° < 9 < 40°; as for the other range 140° < 9 < 165°, these two branching ratios are: 

Br{B- fo{980)K*-) = (7.5 ~ 13.5) x 10^^ (52) 
Br{B^ -> /o(980)i?*°) = (6.7 ~ 12.5) x 10"^ (53) 

where only the central values of other input parameters are used. It is easy to see that the 
pQCD predictions can account for the measured value or the experimental upper limit in 
the range 140° < 9 < 165° (shown in Fig. 2). From the Fig.2(b), one can find the branching 
ratio of B^ — t- /o(980)^*° should be not far away from the upper limit (i.e. 8.6 x 10^^). 
If we take 9 = 140°, the value of Br{B° /o(980)i?*°) is about 6.8 x 10"^ which is 
consistent with the experimental value, (5.2 ± 2.2) x 10~^ But for 25° < ^ < 40°, 

the predicted B^ — )• /o(980)-^*'^ rate exceeds the current experimental limit. 
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FIG. 2: The 9 dependence of the branching ratios (in units of 10^^) of the decays (a) B~ — )• 
/o(980)ii'*~ and (b) — )• /o(980)Er*''. The horizontal sohd hnes show (a) the measured value 
and (b) the experimental upper limit, respectively. The horizontal band within the doted lines 
shows the experimentally allowed region of decay B~ — )• /o(980)i^*^. The vertical bands show 
two possible ranges of 9: 25° < 9 < 40° and 140° <9 < 165°. 



TABLE II: Decay amplitudes for B^ /o(980)i^*" (xlQ-^GeV^), where "this work" denotes 
the results using the distribution amplitudes ^fg, and <I>J^ given in the previous section, 



13l|" denotes the results using the DAs proposed in [l^]. 



ss 










This work 


6.02 


1.12 + 4.37i 


-0.45 - 0.78i 


0.32 + 7.32i 


[13] 


3.25 


0.29 + 0.31i 


0.56 - 0.70i 


-7.49 + 0.42i 


fin 








e/o 


This work 


-6.74 


-19.47 - 59.17i 


2.9 + 11.3i 


0.81 - 0.56i 


L13] 


10.5 


7.54 + 6.97i 


0.27 + 0.28i 


-0.37 + 1.8H 


nn 




^afo 






This work 


0.17 + 0.13i 


0.35 - 6.77i 






m 


0.14-0.07i 


-7.42 + 0.19i 







Our results are larger than the previous pQCD results 131] • Part of the reason is in 
taking the different parameters, for example the decay constant of /o(980). The main 
reason is that the author in [l3| neglected the twist-2 contribution but only used the 
twist-3 distribution amplitude 4>f{x), which is symmetry for x — )■ 1 — x. Taking these 
shapes of distribution amplitude would make the contributions from the non-factorizable 
diagrams (c) and (d) cancel with each other. But here we include the twist-2 distribution 
amplitude and use the asymptotic form of the twist-3 distribution amplitude. In this 
case, the contributions from /q emission non-factorizable diagrams are large. In order to 
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FIG. 3: The 9 dependence of the direct CP asymmetry (in units of percent) of the decays (a) 
/o(980)K*- and (b) /o(980)i^*°. The vertical bands show possible ranges of 6: 

25° <e < 40° and 140° <e < 165°. 



show this character, we Ust the numerical results for different topology diagrams of B — t- 
/o(980)i^*° in Table II. In the table, -^g^^^j^^.o and ^f"^)^^ *o denote as the contributions from 
/o emission (annihilation) factorizable contributions and non-factorizable contributions 
from penguin operators respectively. Similarly, F^*°j:^ and M^*°j^ are the K*'^ emission 
(annihilation) factorizable contributions and non-factorizable contributions from penguin 
operators, respectively. ^e(a)_ft:*o denote the /o emission non-factorizable contribution 
from tree operator O2. It is easy to see that ^f°(j)^*o ^^^d ^/"^^.o obtain an enhancement 
compared to previous estimates. It suggests that the non-factorizable type amplitude is 
sensitive to the shape of the distribution amphtudes. 

Now we turn to the evaluations of the direct CP-violating asymmetries of B~ 
fo{980)K*- and 5° ^ /o(980)ir*° decays in the pQCD approach. The direct CP-violating 
asymmetry can be defined as 

\M\^ + \M\'- ^ 

For the decay — ?► /o(ss)-S'*°, there is no tree contribution at the leading order, so 
the CP asymmetry is naturally zero. But the CP asymmetry of B^ — t- fQ{nn)K*^ is 
large, for the /o emission non-factorizable diagrams (Fig. 1(c) and (d)) give the large tree 
contributions, and its the direct CP asymmetry is about —39%. It is similar to the decay 
B^ — 7- /o(980)i^'*~. From the Fig.3(a), one can find that if taking the mixing angle 
25° < < 40°, the direct CP asymmetry of the decay B' fo{980)K*- is: 

A'^SUB' ^ /o(980)i^*-) = (-15 ~ -25)%, (55) 

which may suffice to explain the experimental result j25| : 

A'^jpiB- fo{980)K*~) = (-34 ± 21)%. (56) 
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But if wc take the mixing angle 140° < 6 < 165°, the value has the opposite sign with 
the experimental result and becomes (10 ~ 33)%. Certainly, the errors from both the 
experimental result and the prediction are large. 

From the upper analysis to the branch ratios of the decays — > /o(980)X*°, B~ — > 
/o(980)ii'*~, it supports the conclusion that the mixing angle should be in the range of 
140° < 6 < 165°. But unfortunately the range of 25° < 6 < 40° as it seems cannot be ruled 
out absolutely. From Fig. 2 and Fig.3, one can find there exist some symmetries for these 
two angle ranges. Within (large) theoretical errors, the results for the two angle ranges 
are both in agreement with the data. For example, if wc take the angle 25° < < 40° in 
the Fig. 3(b), the direct CP asymmetry of the decay B~ /o(980)i^*~ is: 

A'^Sp{B- foi980)K*-) = (0.8 ~ 0.9)%, (57) 

and A'^SpiB' ^ /o(980)fs:*-) = (-1.2 ~ -5.9)% for 140° <9 < 165°. That is to say the 
values of Afi'p{B' ^ /o(980)A'*-) for these two 6 angle ranges are close and both small. 

V. CONCLUSION 

In this paper, we calculate the branching ratios and CP-violating asymmetries of 
-> /o(980)-Sr*° and B~ — )■ fo{980)K*~ decays in the pQCD factorization approach 
by identifying /o(980) as the composition of ss and nn = {uu + dd)/\/2. Using the decay 
constants and light-cone distribution amplitude derived from the QCD sum-rule method, 
we find that: 

• After including the twist-2 distribution amplitude and using the asymptotic form 
of twist-3 distribution amplitude, our results are larger than the previous pQCD 
predictions and can explain the present experimental data or the upper limit. 

• From the results, it indicates that the contributions from the non-factorizable /o 
emission type diagrams are large, at the same time this type amplitude is sensitive 
to the shape of the distribution amplitudes. 

• The branching ratio of S — >■ /o(980)i^* depends on the mixing angle 9 of strange 
and nonstrange components of the /o(980). One can find that there exit some 
symmetries for the values in the two angle ranges (i.e., 25° < 6 < 40° and 140° < 
9 < 165°). So it is difficult to confirm the value of the mixing angle, unless we can 
get enough and precise experimental data. 

• For the neutral decay B^ — )■ /o(980)A'*°, wc predict that the direct CP-violating 
asymmetry is small, only a few percent, which can be tested by the future B factory 
experiments. 
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